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Abstract 

The soft photon theorem in [/(I) gauge theories with only massless charged particles 
has recently been shown to be the Ward identity of an infinite-dimensional asymptotic 
symmetry group. This symmetry group is comprised of gauge transformations which 
approach angle-dependent constants at null infinity. In this paper, we extend the 
analysis to all U{1) theories, including those with massive charged particles such as 
QED. 
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1 Introduction 

The soft photon theorem [1-5] has played a ubiquitous role in the study of QED and more 
general abelian gauge theories. For example, it is essential for taming otherwise uncontrol¬ 
lable infrared divergences in the iS-matrix and is central to the analysis of jet substructure. 
Recent considerations [6-14] have demonstrated that, in abelian gauge theories with only 
massless charged particles, the soft theorem is a Ward identity of an infinite-dimensional 
symmetry group comprised of certain ‘large’ gauge transformations which do not die off at 
infinity. These symmetries are spontaneously broken and the soft photons are the Goldstone 
bosons. This is but one instance of a recently-discovered universal triangle connecting soft 
theorems, symmetries and memory in gauge and gravitational theories [6-10,15-34]. 

Of course in the real world QED has massive, not massless, charged particles. Hence, it 
is desirable to extend our results to the massive case. That goal is achieved in this paper. 
As seen below, the massive case is rather more subtle than the massless one and requires a 
careful analysis of timelike infinity. 

We hope that the identification given herein of the symmetry which controls the electro¬ 
magnetic soft behavior of QED and more generally, the Standard Model will have practical 
utility for organizing and predicting a variety of physical phenomena. 
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The outline of the paper is as follows. In section 2, we establish conventions and review 
relevant aspects of abelian gauge theories and their asymptotic symmetries. In section 3, we 
discuss the asymptotic states, derive the Ward identity of the asymptotic symmetries, and 
demonstrate its equivalence to the soft photon theorem. 

A key ingredient of our analysis is that, in physical applications, the electromagnetic 
field is generically^ not smooth near spatial inhnity Rather it obeys a matching condition 
near which identihes its value at the future of past null inhnity with its value at the 
antipodal point on the sphere at the past of future null inhnity (Xj^). In the appendix, we 
show in detail how this follows from the standard Lienard-Wiechert formulae. 

After completion of this work, we received the eprint [35] by Campiglia and Laddha 
who independently arrive at the same conclusions. They use an elegant method involving a 
natural hyperbolic slicing of Minkowski space. 

2 Abelian gauge theory with massive matter 

We consider the theory of an abelian gauge held coupled to massive matter helds Tj 
with charges eQi, where Qi is an integer, in Minkowski space. In retarded coordinates, the 
Minkowski metric reads 

ds^ = —dt^ + {dx^Y = —du^ — 2dudr + 2r'^jzzdzdz, (2.1) 

where u is retarded time and 7^2 is the round metric on the unit radius 3“^ with covariant 
derivative D^. The S'^ coordinates {z, z) are related to standard Cartesian coordinates by 

r = ZiZ", u = t — r, x^ = rz\z,z). ( 2 . 2 ) 

In retarded coordinates, future null inhnity (X’*') is the null hypersurface (r = 00 , u, z, z). 
Near past null inhnity (X“), we work in advanced coordinates (n, r, z, ^) with line element 

+ 2dvdr + 2r^'^zzdzdz. (2.3) 

Advanced coordinates are given by 

r = Xiz\ v = t + r, x^ = —rxYz,z), (2.4) 

^For instance when, as in electron-positron scattering, the dipole moment is not constant in the far past 
or future. 
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and X~ corresponds to the null hypersurface (r = oo,n, z,f). Note in particular that the 
angular coordinates on X"*" are antipodally related to those on X~ so that a light ray passing 
through the interior of Minkowski space reaches the same value of z at both X”*" and X~. 
We denote the future (past) boundary of X+ by X\ (Xj^), and the future (past) boundary of 
X- by X- (XI). 

We consider theories with a U{1) gauge field strength F = dA subject to the Maxwell 
equation 

= eV,, (2.5) 

where Ji, is the matter charge current. This is invariant under the gauge transformations 

Af,{x) Af,{x) + d^e{x), Ti(a;) i{x), (2.6) 

where e ~ e + 27r and Tj is a wavefunction or field. Gauge transformations that vanish at 
inhnity correspond to redundant descriptions of the same physical state and can be eliminated 
by a choice of gauge. However, as in the massless case [6], we are interested in certain angle- 
dependent large gauge transformations which act non-trivially on physical states. 


2.1 Asymptotics 

We now analyze the behavior of the theory near X'^ in retarded radial gauge 

Ar = 0, = 0. (2.7) 


This gauge choice leaves unfixed a class of residual large gauge transformations parameterized 
by an arbitrary function e~^(z, z) on S'^. These gauge transformations change boundary data 
at X~^ and are to be regarded as physical symmetries of the theory. Near X~^, we assume the 
asymptotic expansion 


n=l 



^‘=11 




u, z, z) 


71=0 


( 2 . 8 ) 


A similar asymptotic expansion holds for fields near X~. 

We are interested in scattering processes for which the initial and final states consist 
of non-interacting massive charges moving at constant velocities. Hence, we require that 
the only contribution to the electric and magnetic fields at future/past timelike inhnity 
±C)o) are those helds sourced by the constant velocity massive charges, and that the 
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magnetic fields vanish at spatial infinity: 


Fzz\x+ = 0) F^ 2 \x- = 0- (2-9) 

In retarded coordinates, Maxwell’s equations read 

r~‘^dr{r‘^Fru) - duF^ + r~‘^{D^F^u + D'^^F^u) = e^Ju, ( 2 . 10 ) 

r-‘^dr{r'^Fru) + r-‘^{D^F,r + D~^F-,r) = ( 2 . 11 ) 

dr{Frz - Fuz) - duFr^ + T-'^D^Fg^ = ( 2 . 12 ) 

Massive particles with hnite energy cannot reach X, so the matter current vanishes at this 
surface: 

Jf.\i = 0. (2.13) 

The leading order equation for the evolution of the gauge held along X^ is then given by 

duF^l^ + du{D^Af + D'^Af) = 0. (2.14) 


( 2 ) 

The free data at this order includes the boundary data FrJ\x+ along with the radiative mode 
Af\u, z, z). 

In advanced coordinates, we can perform the analogous large-r expansion near X“ and 
obtain the leading order equation 

- d,{D^Af + D~^Af) = 0. (2.15) 

The free data at this boundary includes the held strength boundary data Frv | j- along with 
the radiative mode A^z\v^ z, z). The residual large gauge symmetry is parameterized by an 
arbitrary function e~{z,z) on S"^. 

2.2 Matching near spatial infinity 

The above discussion treats the asymptotic dynamics at X+ and X" separately. However, 
to study the semiclassical scattering problem, we must hrst specify how to relate free data 
and symmetry transformations at to their counterparts at X“. Generic solutions to the 
sourced Maxwell equations satisfy^ 

^See appendix for an expanded discussion of this matching condition. 
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F>i>(z,z)\it^F’SHz.z)\j-. (2,16) 

Recalling that, according to (2.2) and (2.4), the points labelled by the same {z, z) in retarded 
and advanced coordinates are antipodally related, this eqnates the bonndary valnes of past 
and fntnre helds at antipodal points near spatial inhnity As discussed in [6,21], a CPT 
and Lorentz-invariant matching condition for the gauge field is given by 


A^(z,z)l^+ = A^{z,z)\ 


XT' 


(2,17) 


Requiring that the large gauge transformations preserve this matching condition gives: 


e*{z,z) = e {z,z). 


(2,18) 


This matching condition singles out a canonical diagonal subgroup of the large gauge trans¬ 
formations at and X~. The corresponding gauge parameters are constant along the null 
generators of X and generate nontrivial physical symmetries of the iS-matrix. 


2.3 Mode Expansions 

The standard mode expansion for the gauge field in the plane wave basis takes the form 

(Pq 1 


A,(u,r,z,z)-eY: j 


(27r)3 2n;/ ^ 






(2,19) 


The free data is contained in the 0{r^) term in this expansion, which we may isolate using 
the saddle point approximation; 




U, Z, Z) = 


le 


POO 

( 2 - 20 ) 

JO 


2 ( 2^)2 

To extract the contribution from the zero modes, we define the following operator: 


z)= du e'-^^duAfiu, z, z) 


22 / )h(a; - cug) -1- e"a„(a;g£)'^5(a; Uq)] . 
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We can separate this operator into its positive and negative frequency components 


a„(a;£), z) =-—^ (2.21) 


Fuz{z,z) = 

a 

with a; > 0 in both expressions. The zero mode is then given by 


FLiz, z) = ^ lim lim ^ [ue*^aa,{ux) + a;£>„(a;x)^] , ( 2 . 22 ) 


Sir 


and creates/annihilates soft photons. An analogous construction holds at X with the in¬ 
coming soft photon operator given by 




uje*'^aa{—(^x) + uefaai—izx) 


(2.23) 


2.4 Lienard-Wiechert Fields 


In the analysis that follows, we will need expressions for the electric field due to moving 
point charges, commonly known as Lienard-Wiechert helds. The radial electric field due to a 
single particle of charge eQ, moving with constant velocity [3 and passing through the origin 
at t = 0 is given by 


Er{t,r,z, z) 


Qe^ 7 (r — tx(z, z) ■ /3) 

4vr _ j.x(^z, z) ■ /S]^ — F + 


(2.24) 


Here x(z, z) is a unit vector specifying a point on the sphere and 7 “^ = 1 — 

The Lienard-Wiechert field near X’*' due to a set of particles, each with charge eQk and 
moving with constant velocity 13k is derived by taking a superposition of the single-particle 
helds (2.24), writing them in retarded coordinates, and taking the large-r limit with u = t — r 
held hxed 


Ft {z, z) = Y^ 


Qke^ 


(2.25) 


fc 47r7^V[l_£(^,7)./5^]2 

Likewise, the held near X“ is derived by taking the large-r limit of the held (2.24) in advanced 
coordinates with hxed v 


Er [z, z) 


Qke^ 1 

V 47 r 7 |r 2 [I + x{z,z) ■ ^kf 


(2.26) 


Importantly, the Lienard-Wiechert formula (2.24) implies that the value of Ey near spatial 
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infinity depends on how it is approached. In particular, and at a fixed angle from 
the origin are not in general equal near rather they obey the antipodal matching condition 
(2.16).3 

For unaccelerated charges, the asymptotic electric held and the asymptotic magnetic held 
B = X X E are time-independent. Since the “hard” radiative photons involved in the scat¬ 
tering process exit/enter at hnite values of retarded/advanced time, the electromagnetic 
helds at and XZ arise solely from the collection of charged particles long after/before the 
scattering process occurs and thus are of the form given above. 

3 Symmetries of the <S-matrix 

In this section we determine the phase associated to a large gauge transformation on a 
(dressed) asymptotic massive charged particle state, hnd the iS-matrix Ward identity and 
hnally demonstrate its equivalence to the soft photon theorem. 

3.1 Gauge transformations of asymptotic states 

Outgoing massless particles of charge eQ and momentum p, as considered in [6], pierce X'^ 
at a dehnite point (z(p),^(p)). The associated out-state therefore acquires a phase 

(3.1) 

under a large gauge transformation. Here we are interested in massive particles that never 
reach X"*", so determining the associated phase is more subtle. There is no canonical point on 
the associated to a massive particle^. Indeed a massive particle with zero three-momentum 
is rotationally invariant. In this subsection, we use the Lienard-Wiechert formula to deter¬ 
mine the analog of the phase (3.1). The asymptotic states associated to the QED iS-matrix 
are typically taken to be free photons and “bare” non-interacting charged particles. However, 
infrared (IR) divergences entering 1-loop calculations render such matrix elements formally 
inhnite. Unlike ultraviolet divergences, which are associated to high energy, unobservable 
physics, these IR divergences have a clear physical interpretation associated to low energy 

^The constant-velocity trajectory considered in this section gives rise to a Lienard-Wiechert field that 
is insensitive to the choice of Green’s function. In the appendix, we consider slightly more complicated 
trajectories to demonstrate that this matching condition holds for a generic Green’s function. 

^Since the wavefunction of a massless particle localizes to a point on the conformal sphere at null infinity, 
it is related to a local operator insertion on that sphere. Massive particles, on the other hand, will involve 
non-local weighted integrals of operators over the sphere. 
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observables. This well-understood problem stems from an improper choice of asymptotic 
scattering states [36-43]. Formally, the standard assumption that the asymptotic dynamics 
is governed by the free Hamiltonian Hq does not hold in QED. It fails precisely because 
the long-wavelength, infrared sector of the gauge theory fails to decouple from the charged 
particles at early/late times. Physically, it is impossible to separate a charged particle from 
its electromagnetic held. The ability to do so would violate Gauss’ law, which can be imple¬ 
mented as a constraint in the quantum theory: 


V • E{x, t) 


e^p{x,t) 


\phys) = 0. 


(3.2) 


The resolution, worked out in [36-43], is to choose scattering states that diagonalize the 
asymptotic Hamiltonian Has, dehned to account for the long-range interactions in QED. 
In particular, “undressed” charged particles without electric helds are not eigenstates of the 
asymptotic Hamiltonian and cannot appear as scattering states. Consider the case of a single 
“undressed” electron state 

V^(x)lO). (3.3) 

This state is not invariant under small gauge transformations and does not satisfy the Gauss 
constraint. To form a physical state, we must also include the electron’s electric held. This 
problem can be solved by introducing dressed states of the general form [36,44-46] 


|V'(x),E(x)) ='0(x)exp — 


Ei{x,y)A\y)d^y 


| 0 ) 


for any Ei{x,y) satisfying the Gauss constraint 


(3,4) 


V‘Ei(x,y) = Qe^S^(x-y). 


(3.5) 


Under a general gauge transformation, this state transforms as 


(x), E{x)) —)■ exp 


'52 


(fzjzzr'^e{z,z)Er{z,z) 


(x),E(x)). 


(3.6) 


For small gauge transformations e{x) vanishing at inhnity, the state is left invariant, while 
large gauge transformations rotate the state by a phase. 

Suppose that the outgoing states with momentum have been chosen to 

diagonalize the asymptotic Hamiltonian and are dressed with their Lienard-Wiechert helds 








so that the electric held at TJt is given by 


Er{r,z,z)\p)out = 


Qe^ 


[1 _ x{z, z) ■ /3]2 

Under a large gauge transformation, such a state acquires a phase 


\P) 


out' 


\p)out -t exp 


i / d 


Q 


/x+ 


47 r 72 [1 _ x{z, z) ■ /3]2 


\P) 


out • 


Similarly, in-states transform as 


\p)in -)> exp 


i / d^z^^2£~ 


Q 


[l+x{z,z) ■ PY 


\p)i 


(3.7) 


(3.8) 


(3.9) 


For an n-particle state, the phase will be a sum of n such terms. This phase replaces the 
much simpler expression (3.1) for massless particles but nevertheless, as will be seen shortly, 
precisely reproduces the soft factor for massive particles. 

In principle, a state could be dressed with a different solution to the Gauss constraint. 
However, any such solution can only differ from the Lienard-Wiechert solution by a solution 
to the source-free Maxwell equations. These more general radiative solutions are eliminated 
by our boundary conditions at and XZ, which single out the sourced, non-radiative 
contribution to the held. Physically, we want to dress our states with a ‘minimal’ photon 
cloud that satishes the Gauss constraint but excludes additional photons corresponding to 
radiative contributions to the electromagnetic held. These radiative contributions cross X 
and do not pass through 


3.2 Ward Identity 

We are now in a position to discuss the symmetries of the iS-matrix. The symmetry trans¬ 
formations (2.6) for massless matter helds have already been analyzed in [6], where it was 
demonstrated that the charge 

Qt = \ [ d‘^z^,^e^{z,z)Fl:l\z,z) (3.10) 

e Jx+ 

generates the correct X”*" symmetry transformation on the gauge held and matter helds. The 
form of this charge is essentially hxed by the transformation law for the gauge held. We 
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can use the leading order Maxwell equation (2.14) to turn this expression into an integral 
over X+. As discussed in section 2.1, the existence of massive particles generates charge flux 
through future timelike inhnity, so the local charge operator takes the form 

Qt = \ f Izzducfz e'^du{D^Af'> + + \ f (fz 

e' Jx+ e2 J^+ 

= 4 / + D^F^,) + 4 / dh lz,e+Fl:ll (3.11) 

Js^ Jx+ 

The hrst piece of the charge is written in terms of the soft photon operator and will 
be referred to as the soft charge Qg. If we consider the fixed-angle charge by choosing 
e{z, z) = d^(z — w), then the second term is simply the radial electric held in the direction 
{w,w) resulting from the charged particles in the hnal state at We label this term the 
hard charge It differs from the expression for the hard charge in the massless case which 
involves an integral over X”*". 

An analogous computation can be performed at X“, where the hard charge encodes the 
radial electric held of the charged particles in the initial state at i~. The charge is given by 

\ f e-(lTF’i + + 4 / = Qs + Qh- ( 312 ) 

Js^ Jxz 

The statement that the transformations (2.6) are symmetries of the iS-matrix is equivalent 
to the statement that the charges (3.11) and (3.12) commute with the iS-matrix: 

{out| — SQ~^ |in) = 0. (3.13) 

In order to facilitate comparison with the soft theorem, we separate the hard and soft 
contributions and rearrange the Ward identity: 

(out| {Qs<S - SQg) |in) = -(out| {QjjS - SQ]^) |in). (3.14) 
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3.3 Soft Theorem Ward Identity 

The soft photon theorem for the emission of an outgoing photon in a scattering process with 
m incoming hard particles and {n — m) outgoing hard particles reads 


lim 

u—>-0 


\aa{q)S\pi,...) 


eco 


n 


Z <3“ 

_k=m-\-l 


Pk * 

Pk • q 


m 




Pk * ^OL 
Pk-q 




(3.15) 

A null momentum vector is uniquely specihed by an energy and a point z on the asymptotic 
sphere, and so we parameterize the photon’s momentum as 


= uj[l,x{z,z)] = uq^{z,z), 


(3.16) 


where £ : —?■ is an embedding of the sphere into flat three-dimensional space. 

We parameterize a massive particle’s momentum as 


Pk = 


(3.17) 


where m is the rest mass of the particle, (3 is the particle’s velocity and 7 is the relativistic 
factor 7“^ = 1-/1^. 

We can relate the left-hand side of equation (3.15) to the zero mode operator defined in 
equation (2.22) by taking a weighted sum over polarizations. If we perform the analogous 
operation on the right-hand side and use the identity 


d:,x\z,z)^i 


Y Pk * 

Pk ■ q{z,z] 


= dJog{pk ■ q), 


(3.18) 


the soft theorem can be written 


(p™+i,...|F°^5|pi,...) = - 


Stt 


Qk dz \og{pk ■ q) - '^Qk log(pfc • q) 


_k=m-\-l 


k=l 


(Pm+l,...|5|pi 

(3.19) 


The soft photon theorem for an incoming soft photon reads 


lim u{pra+i, ■ ■ ■ |5a„(g)^|pi,...) = 

ui^O 


—eu) 




fc=m+l 


Pk ■ q 


k=l 


Pk-el 

Pk ■ q 


{P-. 


m+l) 


|5|pi,. . 

(3.20) 
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An identical calculation yields 


{Pm+l,. . .\SF^^\pi,. . .) 


Stt 


n m 

Qk dz log{pk ■ q') Qk \og{pk ■ q') 

_k=m-\-l k=l 




( 3 . 21 ) 


where q' = [1, —F{z,z)]. Taking the divergence of each equation, using global charge con¬ 
servation, and integrating against the respective gauge parameter, we hnd 


(p™+i. 5|p„ ...) 

= -i [£'+Ji„)(p„+i....|5|pi,...) (3.22) 

and 

(p„+i,... |5 U iez^.,e-{D‘FS, + d^fS,)] Ipi, ...) 

= i j- [Brim)(Pm+i,... |S|pi,...). (3.23) 

Taking the difference and using the matching conditions (2.16)-(2.18), we obtain 

{Pm+i,---\Qs^-<SQs\pi,...) 

=-^ j (fz'y;,sr‘^(^e^Er\^+-e~Er\^.'^{pm+i,---\S\pi,...) (3.24) 

= ~{Pm+i, ■ ■ ■ \Q~h^ ~ ^Qh\p^^ ■ ■ ■)■ (3.25) 

This precisely reproduces the Ward identity (3.14). 

In conclusion, while the details are more intricate than the massless case, the soft photon 
theorem is the Ward identity of an inhnite-dimensional asymptotic symmetry group for 
abelian gauge theories with massive particles. We expect similar conclusions apply to other 
contexts such as non-abelian gauge theory and gravity. 
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4 Appendix: Gauge field strength near 


In this section, we consider an idealized semiclassical scattering process in which m incoming 
massive particles with constant velocities {/3i,... ,/3m} scatter to (n — m) outgoing massive 
particles with constant velocities {fim+ii ■ ■ ■ For scattering occurring at the origin at 

t = 0, the semiclassical electromagnetic current is given by 

m p Up 

Ux) ^Y,Qk dTU^e(-T)6^‘‘>(x - U^t) + 5^ qO dTl7‘e(T)i<‘)(i - f/V), (4.1) 

k=l k=m+l 


where = 7fc[l, Pk] is the 4-velocity of the particle. Ignoring the radiative contributions 
arising from the inhnite acceleration of particles at the origin, the held strength sourced by 
this current takes the form 


Frt{x) 


^ g{x; Pk) Qklkjr - tx • /3fc) ^ -A h(x; jSk) Qkjkjr - tx • /3fc) 

k=i l7|[A - rx ■ /3fc]2 - F + r2|3/2 47r \-fl[t - rx ■ ^kf - F + ’ 

(4.2) 


where the functional form of g and h depends on the choice of Green’s function. 
For the retarded solution, the asymptotic behavior of g and h is given by 


5r(r = oo,u,x;/5fc) = 0(-u), h(r = oo, u, x;/3fc) = 0(u), (4.3) 

^(r = cx),u,x;/3fc) = 1, h(r = cx), u, x; ^a,) = 0. (4.4) 


The electric held at is obtained by working in retarded coordinates and taking the limit 
r —)■ oo, followed by the limit u —?■ —oo. This electric held will be of the form (2.25), but 
only receives contributions from the incoming particles. On the other hand, the electric held 
at X^ is obtained by working in advanced coordinates and taking the limit r ^ oo, followed 
by the limit u —)■ +cxo. The electric held measured at will be of the form (2.26), but will 
also only receive contributions from the incoming particles, thereby satisfying the matching 
condition (2.16). 
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Likewise, for the advanced solution, the asymptotic behavior of g and h is given by 

g{r = oo, u, X] /3k) = 0, h{r = oo, u, x; /3k) = 1, (4.5) 

g{r = oo,v,x; /3k) = Q{-v), /i(r = oo, n, x;/3fc) = 0(n). (4.6) 

Hence, the advanced solution also obeys the matching condition (2.16) near but in contrast 
to the retarded solution, only receives contributions from the outgoing particles. Moreover, 
linear combinations of the advanced and retarded solutions evaluated near will obey the 
matching condition and receive contributions from both outgoing and incoming particles. 

Of course, one could always add a homogeneous solution to the free Maxwell equations 
which does not obey the matching condition. However, we do not know of any physical 
application in which it is natural to do so: hnite energy wave packets die off at Hence, we 
conclude that the antipodal matching condition (2.16) holds in generic physical applications. 
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